Crossover between Hydrodynamic and Kinetic Modes in Binary Liquid Alloys 
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Inelastic x-ray scattering (IXS) measurements of the dynamic structure factor in liquid ^57X43, 
sensitive to the atomic-scale coarse graining, reveal a sound velocity value exceeding the long wave- 
length, continuum value and indicate the coexistence of two phonon-like modes. Applying Gener- 
alized Collective Mode (GCM) analysis scheme, we show that the positive dispersion of the sound 
velocity occurs in a wavelength region below the crossover from hydrodynamic to atom-type ex- 
citations and, therefore, it can not be explained as sound propagation over the light specie (Na) 
network. The present result experimentally proves the existence of positive dispersion in a binary 
mixture due to a relaxation process, as opposed to fast sound phenomena. 

PACS numbers: 67.40.Fd, 63.50.+X, 67.55. Jd, 61.10.Eq 



I. INTRODUCTION 

Alkali metals in their liquid state are outstanding ex- 
amples of simple liquids, as they encompass most of the 
physical properties of complex fluids without system spe- 
cific complications. For this reason, in the last fifty years 
they have been the subject of intensive investigations, 
both experimentally and by means of computer simula- 
tions, aiming at an understanding of the mechanisms un- 
derlying their collective dynamics at a microscopic length 
scale 1 * 2 -. Interest in the collective excitations and inter- 
diffusion processes in liquid binary mixtures, was stimu- 
lated by the pioneering Molecular Dynamics (MD) study 
on the dynamics of the liquid Na-K alloy by Jacucci 
and Mc Donald 3 -. Disparate masses metallic alloys soon 
became the subject of intensive studies: the existence 
of a new, high frequency mode with phase velocity ex- 
ceeding the hydrodynamic value was reported by MD 
simulations^,^, and experimentally proved by Inelastic 
Neutron Scattering (INS) in Li 4 Pb ; 7 -^, in NaCsi^ and 
Li3oBi7o^i. The very nature of this additional mode, tra- 
ditionally named fast sound, has been in focus of a lively 
debate in the last three decades. It was either inter- 
preted as an atom-type acoustic mode, supported by the 
light (Li) ions only, which would merge at low Q with 
the corresponding low frequency mode (Pb) into a sin- 
gle hydrodynamic velocity, or like an optic-like excitation 
arising above a certain treshold Q value. The existence 
of a crossover from hydrodynamic into atom-type excita- 
tions can be theoretically rationalized within a concept of 
kinetic (non-hydrodynamic) modes in the framework of 
Generalized Collective Modes (GCM) approac h 12 i 13 . It 
has been shown that, at low Q, two modes exists and can 
be associated to hydrodynamics density and concentra- 
tion fluctuations (collective region) while, above a certain 
characteristic wavenumber (dependent on mass ratio), 
each of the two excitation reflects the dynamics of a single 
atomic specie (atom- type excitations). The experimental 
identification of the crossover, however, has been heavily 
debated. This was mostly due to the inherent difficulty 
of disentangling the coherent contribution from the to- 



tal INS scattering signal, paired by the kinematic limita- 
tion restricting the accessible energy-momentum region. 
In He-Ne mixtures, for instance, the upper limit of the 
hydrodynamic region has been questioned, also in view 
of the possible definitions of the excitation frequency in 
terms of the maxima of the Dynamic Structure Factor 
(DSF), S(Q, ui), rather that those of the longitudinal cur- 
rent autocorrelation functio n 14 i 15 i 16 ' 17 . 

C l (Q,u;) = ^S(Q,lu) (1) 

Similar controversies appeared in the case of 
watei— — — — — . where for long time the existence 
of fast sound was debated, until the viscoelastic origin 
of such phenomenon was clarifie d 23 ' 24 ' 25 . In fact, 
the advent of the new radiation sources opened the 
possibility to perform Inelastic Scattering with X-rays 
(IXS), overcoming the previously mentioned limitations 
of INS. The purely coherent X-ray atomic cross section 
on one hand, and the relative high energy of the prob- 
ing X-rays on the other, allowed to study the purely 
coherent dynamics of disordered systems over a wide 
energy-momentum range covering the region across the 
first pseudo Brillouin region (i.e. momentum transfers 
up to the first maximum of the static structure factor, 
Qm)- As a consequence, the last decade saw a renewed 
interest in the investigation of high frequency dynamics 
in liquids and glasses. 

It is worth to emphasize in this context the lack of 
a unique and reliable methodology to analyze scattering 
experiments performed on binary liquids. In majority of 
experimental studies one makes use of either single DHO 
model 2 ^ or even a memory-function ansatz -designed for 
one-component liquids- to estimate the dispersion law 
of collective excitations in binary melts2£. Within this 
kind of approaches one necessarily ignores the existence 
of non-acoustic high-frequency excitations in binary liq- 
uids, which beyond hydrodynamic region contribute to 
the shape of partial dynamical structure factors. In this 
study we show how one can apply more sophisticated 
theoretical schemes to the analysis of experimental data 
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in a molten metallic alloy, namely a theoretical GCM 
approach, which consistently treats contributions from 
hydrodynamic and kinetic modes to dynamical processes 
in liquids. 

We report here on an IXS study of the microscopic dy- 
namics in the molten alkali alloy Na^Ki?, which allowed 
us to: i) Identify the presence of two phonon-like modes, 
an high frequency and a low frequency one. ii) Point out, 
by applying GCM theory to a binary mixture, the origin 
of these two modes and how they contribute to the mea- 
sured IXS spectra, i.e. essentially to the mass density 
fluctuations spectra. Specifically, we show the existence 
of a crossover from a hydrodynamic regime, where only 
one of the two modes is active, to an atom-type regime, 
where two excitations appear in the S(Q,u>). Hi) Assign 
the sound velocity excess over the hydrodynamic value to 
a relaxation process similar to that observed in simple 2 , 
molecular 28 and hydrogen bonding liquids—, as opposed 
to the idea of acoustic excitation propagating over the 
network supported by the light specie or to the effect of 
optic-like excitations. 

II. THE EXPERIMENT 

The experiment reported in this work was carried 
out at the high resolution beam line ID16 of the Euro- 
pean Synchrotron Radiation Facility (Grenoble, Fr). The 
backscattering monochromator and analyzer crystals, op- 
crating at the (11,11,11) silicon reflections gave a total 
energy resolution of 1.5 meV, while energy scans were 
performed by varying the temperature of the monochro- 
mator with respect to that of the analyzer crystals. A 
five analyzers bench allowed us to collect simultaneously 
spectra at five different values of constant momentum 
transfer, covering a Q region below the position of the 
main diffraction peak {Qm «18nm _1 ). The sample con- 
sisted of an NaK alloy at T — 300 K with sodium concen- 
tration Cn s = 57 atomic %, and was hold in a silica cap- 
illary of 2 mm inner diameter and 10 /im wall thickness, 
sealed and kept in an inert atmosphere. Energy scans 
in the range —40 < E < 40 meV where repeated up to 
a total integration time of 300 s/point. Energy scans of 
the empty cell were collected in the same Q — E range 
as for the sample, prior to filling the capillary, and were 
subtracted to the final IXS spectra. The signal measured 
in an IXS experiment is related to the double differential 
cross-section depending on the exchanged momentum Q 
and energy E — huj, and it is basically the convolution 
of the instrumental resolution, R(u>), with the classical 
X-ray weighted DSF S IX s(Q, w), 

I{Q,w) = J <ko'R{u - u')S IXS [Q, a/) % ^/*% T 

(2) 

in which the last term accounts for the detailed balance 
condition. The outcome of the IXS experiment (open 
circles) is reported in Fig. Q] for selected values of the 
exchanged momentum Q. 
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FIG. 1: Selection of IXS spectra from the NaK liquid alloy 
(cjva = 57at%, T=300K) measured for fixed Q values (nm -1 ). 
The continuous line on the anti-Stokes side displays the out- 
come of GCM theory, while on the Stokes side shows Eq(6]in 
which the eigenmodes are obtained as best fit to the data, as 
explained in the text 



In a binary mixture, Sixs(Qj^) can be represented 
as a linear combination of the partial dynamic structure 
factors 



S ij (Q,u) = J dte~ iut (m(Q,0)rij(Q,t)) (i,j = Na, K) 

(3) 

with weights depending on the Na and K atomic form 
factors as well as on the concentration 2 ^. In the pre- 
vious relation rii(Q,t) is the number density of species 
i. Alternative representations of Sixs(Qt^) can be 
given in terms of dynamical variables derived from the 
n.i as independent linear combinations. The most used 
options, in this respect, are the Bathia Thornton 2 - 
number-concentration, or the March mass-concentration 
variable^. 
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III. THEORETICAL MODEL 

The usual approach to describe the DSF in a liquid 
is based on the solution of the Langevin equation for 
an appropriate set of dynamical variables. This can be 
done, for instance, by modelling a suitable memory ker- 
nel. At least two relaxation times must be accounted 
for in the second order memory function to describe the 
density fluctuations in a monoatomic liquid 2 ", and a non- 
trivial generalization for the case of binary liquids is re- 
quired, along the line of the scheme proposed in Ref£. 
For this reason we have chosen the GCM approach, which 
additionally treats coupling with thermal fluctuations 
and permits to solve the generalized Langevin equation 
(GLE) in terms of dynamical eigenmodes^ 3 -. 

In such an approach the GLE is solved in a markovian 
approximation, while correlations faster than the hydro- 
dynamic ones are accounted for by an appropriate choice 
for the set of dynamical variables. On this ground the 
solution of the GLE is straightforwardly given in terms of 
eigenvalues and eigenvectors of the hydrodynamic matrix 



T(Q)=F(Q,t = 0). F(Q 



0) 



(4) 



where F(Q, t) is the matrix of time correlation functions 
of all the variables belonging to the set, and F(Q,z) is its 
Laplace transfor m 12 ! 13 . The partial DSFs Sij(Q,oj) are 
then obtained by means of a time Fourier transform of 
the corresponding components F i j{Q 1 t) related, for ex- 
ample, to the partial number densities m and rij. As we 
will show, the GCM approach captures all the essential 
features of the reported IXS experiment. 

A suitable extended set of hydrodynamic variables for 
a binary mixture is provided by the se%^ 

A (8) = {n t (Q, t),n x (Q, t), J t (Q, t), J X {Q, t),s(Q, t), j t (Q, t) 

(5) 

where the hydrodynamic variables reflecting the slow- 
est fluctuations in the binary liquid are: total number 
density nt(Q, t), mass concentration density n x {Q, t), to- 
tal longitudinal mass current Jt{Q, t) and energy density 
e(Q,t). The variables with overdots correspond to the 
first time derivatives of the relevant fluctuations and are 
introduced to give a better representation of those pro- 
cess that are faster than the hydrodynamic ones. It is 
worth to emphasize here that, since total density and 
mass-concentration density can be easily represented as 
linear combinations of partial densities, one can represent 
the set via partial dynamical variables. Hence, the 
generalized hydrodynamic matrix T(Q) constructed on 
t-x dynamical variables © or on partial dynamical vari- 
ables will have identical eigenvalues, because both sets of 
dynamical variables are connected by linear transforma- 
tion. 

Within this scheme, the IXS weighted DSF can be ex- 
pressed as: 



Sixs(Q,u) 



-Re 



E 



GL XS (Q) 
z + z a (Q) 



Both the eigenvalues (z a (Q), which are complex or real in 
the case of propagating or diffusive modes, respectively) 
and the corresponding eigenvectors (G I a xs (Q) which are 
the X-ray weighted eigenvectors, i.e. are determined as 
the appropriate linear combination of partial Na and K 
densities, accounting for form factors and atomic concen- 
tration), can be evaluated within GCM theory once the 
Q-dependent elements of the hydrodynamic matrix are 
known. This requires the knowledge of initial values of 
time correlation functions of the variables set, as well as 
the correlation times of the kind 
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F mn (Q,0) 



F mn (Q,t)dt 



(7) 



with m,n = nt,n x , e—. To evaluate these latter we 
performed molecular dynamics simulations on a sys- 
tem of 8000 particles interacting via effective two-body 
potentials^ in a microcanonical ensemble at T = 298K. 
The production run was performed over 3x 10 5 timesteps, 
and all the static and time correlation functions needed 
for GCM analysis were directly calculated in MD. The 



minimal wavcnumber reached was 0.81 



The re- 



(6) 



suiting spectrum consists of eight dynamical eigenmodes: 
complex-conjugated pairs of eigenvalues correspond to 
phonon-like collective excitations z a (Q) = <r a (Q) ± 
iu/ a (Q) with u> a (Q) and cr a (Q) being frequency and 
damping of a-th excitation, respectively, while real eigen- 
values represent purely diffusive relaxation processes. 

The prediction of parameter-free GCM theory (red 
lines) are reported in Fig. [1] along with the lineshape 
of Eq[5] in which the eigenmodes are obtained through a 
best fit to the experimental data (blue lines). In this lat- 
ter case we used as initial values those from GCM, and 
kept fixed all the frequencies of the propagating eigen- 
J<dMv&) (&t@g&)}iry part of the complex eigenvalues), and 
the (real) eigenvalues of the diffusive modes (diffusion 
coefficient). The weight coefficients G I a xs (Q) as well as 
the damping coefficients (real parts) of propagating com- 
plex eigenvalues were instead let as free parameters. Sum 
rules where used as constraints to further reduce the 
number of free parameters^. The stage of fitting pro- 
cedure is necessary since the GCM analysis is based on 
MD simulations, which use effective two-body potentials. 
Namely the use of effective potentials can introduce some 
extra error, because it is well known, that effective two- 
body potentials very rarely can yield the correct melting 
point, therefore under- or overestimating the contribu- 
tions from relaxation processes to the shape of dynam- 
ical structure factors. In the anti-Stokes side of Fig. [1] 
one can see, that the parameter-free GCM approach cor- 
rectly predicts the frequency of collective excitations ob- 
served in the scattering experiments, only relative weight 
of relaxation processes (central peak) and propagating 
modes is not well reproduced. Namely at that point the 
aforementioned fitting procedure was applied, and later 
we will show, that the fitting procedure does not alter 
the main results of GCM analysis. Figure shows the 
imaginary part of the two complex eigenvalues of the 8- 



FIG. 2: Experimental dispersion relation obtained from the 
maxima of Cixs(Q,u) (x), along with the hydrodynamic 

dispersion derived by ultrasound measurements 32 ( ). 

The two modes predicted by complete 8-variables GCM the- 
ory (high frequency o, low frequency •), and the outcome of 
partial 3-variables GCM analysis with total (blue ) con- 
centration (pink - - -), Na (red — — ) and K (green - — ) 
variables are also reported. 

variables GCM treatment (circles and dots). In order 
to ascertain the nature of these phonon-like modes, we 
performed additional GCM analysis using four projected- 
out onto processes of different origin 3- variable eigenvalue 
problems using basis sets 

A (3t) = {n t (Q,t),J t (Q,t),j t (Q,t)} 

A {3x} = {n x (Q,t),J x (Q,t)J x (Q,t)} 

A^ 3 ^ = {n Na {Q,t) : J Na (Q,t),j Na (Q,t)} 

A (3K) = {n K {Q,t),J K (Q,t),J K (Q,t)} (8) 

This approach allows to ascertain the origin of each 
branch of collective excitations in the spectrum and pro- 
cesses responsible for them in different Q-region a 12 i 33 ' 34 . 
The result is reported in the same Fig. [3] and clearly 
shows how the atom type excitations (dotted blue and 
pink lines for the Na and K subsets, respectively) corre- 
spond to the two eigenvalues of the 8- variable treatment 
in the high Q limit, while at low Q the same eigenval- 
ues are reproduced by the total density and concentra- 
tion subsets. This is a clear indication of the existence 
of two dynamical regimes, a low Q, collective region, 
and an high Q, atom-type region, in agreement with re- 
cent observation of optic-like modes above ~ 2 nm" 1 
in MD simulations on Li alloys 5 . Interestingly, the ex- 
perimental sound velocity q(Q) — toi(Q)/Q determined 
by the maximum u>i of the raw current autocorrelation 
function l) 2 Sixs(Q,u) (black crosses) turns out to be 
mainly determined by the low energy mode at low Q 
and by the high energy mode at high Q. In addition 
to that, it does not reach, even at the lowest explored 
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FIG. 3: Classical resolution deconvoluted S(Q,ui) (continu- 
ous line) and C l {Q,uj) (dashed lines) obtained from the IXS 
measurement. The presence of a single mode exceeding the 
adiabatic frequency (marked by the arrow) is clearly visible 
at Q=2 nm -1 , i.e. well below the crossover to atom-type dy- 
namics. For Q=6.5 nm -1 two excitations appears, with the 
high frequency one dominating at Q=9.0 nm -1 . 

which we report the DSF, S(Q,uj) and longitudinal cur- 
rent spectra, C l (Q,lu), from Eq.© adjusting the weight 
of the high frequency and low frequency modes to the 
IXS spectra. At the lowest Q values a single mode domi- 
nates the spectrum, with an energy clearly exceeding the 
one expected from hydrodynamic dispersion u — cqQ. 
Around the first crossover (Q=6.5 nm" 1 ) both the high 
frequency and the low frequency mode do contribute to 
the spectra, while at Q = 9 nmT 1 most of the IXS signal 
is due to the high frequency, Na-like excitation. 

This picture is substantiated by direct inspection of the 
MD data. Different time correlation functions derived in 
MD simulations and analyzed by the GCM approach can 
reveal different types of collective excitations existing in 
binary liquids. In Figs. [HandElwe show MD-derived lon- 
gitudinal total and mass-concentration autocorrelation 
function at three lowest wavenumbers sampled in sim- 
ulations. It is clearly seen, that the mass-concentration 
current autocorrelation functions C xx (Q,t) functions re- 
flect propagating excitations with much smaller time 
scale, than the sound excitations, which as in the case 
of single-component liquids, contribute to the shape of 
Ctf.(Q,t) functions. Furthermore, almost identical shape 
of C xx (Q,t) functions at small wavenumbers implies, 
that the short-time excitations have almost the same fre- 
quency and damping in the long- wavelength limit, that is 
usually the evidence of kinetic optic-like excitations The 
presence of the two phonon-like modes in the IXS spec- 
tra, at intermediate Q values, can be conveniently quan- 
tified looking at the relative weights reported in Fig[|5J 
The low frequency mode is clearly dominant below the 
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FIG. 4: Normalized total mass current autocorrelation func- 
tion at the three lowest Q values investigated in MD 
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FIG. 5: Normalized mass-concentration current autocorrela- 
tion function at the three lowest Q values investigated in MD 



sharp crossover occurring at Q w 6 nm . Around this 
value both modes contribute to the IXS spectra, while 
above the high frequency mode dominates up to a new 
crossover at Q « 13 nm -1 . The IXS cross section is 
indeed roughly proportional to the total density auto- 
correlation function, and hence below the crossover it 
reflects the collective longitudinal excitation and not the 
optic-like mode. As previously mentioned, the ampli- 
tude values from GCM and from fitting procedure are 
quite different, but the oscillating behavior with the pres- 
ence of sharp crossovers is captured by both the methods. 
The fact that positive dispersion observed from the IXS 
data occurs well below the first crossover, when the spec- 
tral frequency is solely reproduced by the low frequency 
mode, allows us to rule out a "fast sound" mechanism 
i.e. the transition from the hydrodynamic regime to a 
dynamical regime in which the sound propagates over a 
network formed by the light component of the mixture, 
as it happens in disparate masses alloys. On the contrary, 
the origin of such an effect has to be traced back to a re- 
laxation process of viscoelastic origin, similar to that in- 
voked in simple- and molecular 28 liquids, glasse s 35 ' 36 and 
in the controversial case of water and hydrogen bond- 



FIG. 6: High (circle) and low (triangle) frequency modes from 
GCM (red) and from fitting procedure (blu). The crossover 
at Q « 6 nm -1 can be observed. 
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Along the same line, the inspec- 
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FIG. 7: MD calculated C^ aNa (Q,oj) ( ) and C£ K (Q,w) 

( ) for Q=2nm _1 . The same acoustic mode turns out to 

dominate in both species, and no sign of fast sound-like phe- 
nomenology, in which the light component is dominated by 
the high frequency mode, is recovered at Q lower than the 
crossover value (Q ~ 5nm _1 , see Fig. [8} from MD. 



tion of the MD derived spectral functions C^ aNa (Q,uj) 
and Ck K (Q,oj) presented in Fig. [7] for a Q = 2 nm _1 - 
i.e. well below the crossover- is particularly enlightening. 
The fingerprint of fast sound would be here the presence 
of an extra high frequency mode in the light component 
(Na) spectra which, in turn, is absent in the correspond- 
ing heavy component (K) spectra. As can be observed 
in Fig. this is certainly not the present case. Con- 
versely, the same acoustic mode dominates in both the 
species, accompanied by a common weak (Q 2 vanishing) 
high frequency feature. As shown in Fig. [5J this scenario 
holds at any Q's below the crossover, while at high Q's 
the partial spectra are dominated by distinct optic-like 
modes which can be identified with the high and the low 
frequency GCM modes. 
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FIG. 8: MD calculated dispersion relation for Na ( ) and 

K ( ), the high (o) and low (•) frequency GCM modes are 

also reported. At Q's below the crossover the same scenario 
depicted in Fig. [7] holds for the partial Cj^ aNa (Q,io) and 
Cj( K (Q,lj) currents, which are dominated by what can be 
addressed to as the GCM low frequency mode. 



liquid binary mixture Na57lQ3 which suggest the pres- 
ence of two distinct phonon-like excitations beside diffu- 
sive modes. A quantitative analysis based on GCM the- 
ory allows to identify the nature of these excitations. The 
relative weight of the phonon-like modes shows an oscilla- 
tory behavior with Q, ruled by sharp crossovers defining 
collective and atom-type dynamic regions, in which the 
global density/concentration fluctuation and the partial 
Na/K density fluctuations dominate, respectively. The 
observed positive dispersion of the acoustic branch, i.e. a 
sound velocity value exceeding the long wavelength limit, 
occurs well within the collective region. This indicates 
for this process an origin due to the existence of a mi- 
croscopic relaxation, ubiquitous in monoatmic liquids, as 
opposed to the fast sound phenomena observed in dis- 
parate masses mixtures, where enhanced sound velocity 
has been suggested to be either due to the propagation of 
sound over the light component network or to the pres- 
ence of an optic branch. 



IV. CONCLUSIONS 
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